Inversion mechanism for the transport current in type-II superconductors 
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The longitudinal transport problem (current is applied parallel to some bias magnetic field) in 
type-II superconductors is analyzed theoretically. Based on analytical results for simplified configu- 
rations and relying on numerical studies for general scenarios, it is shown that a remarkable inversion 
of the current flow in a surface layer may be predicted under a wide set of experimental conditions. 
Strongly inhomogeneous current density profiles, characterized by enhanced transport towards the 
center and reduced or even negative values at the periphery of the conductor are expected when the 
physical mechanisms of flux depinning and consumption (via line cutting) are recalled. A number 
of striking collateral effects such as local and global paramagnetic behavior are predicted. Our geo- 
metrical description of the macroscopic material laws allows a pictorial interpretation of the physical 
phenomena underlying the transport backfiow. 

PACS numbers: 74.25.Sv, 74.25.Ha, 41.20.Gz, 02.30.Xx 



I. INTRODUCTION 

Typc-II superconductors under the action of a trans- 
port current and a longitudinal magnetic field may ex- 
hibit the counter intuitive phenomenon of negative re- 
sistance within a certain set of experimental conditions. 
This property, together with other intriguing phenom- 
ena, such as the observation of paramagnetic moments, 
and the compression of the transport current by the 
action of a parallel magnetic field, have been reported 
in the course of intense experimental and theoretical 
activities ji"— Most of these works were primarily con- 
cerned with the arrangement of the macroscopic current 
density J along the so-called nearly force free trajecto- 
ries. Recall that if J is nearly parallel to the magnetic 
induction B, moderate or weak pinning forces are needed 
for avoiding the detrimental flux-flow losses related to 
the drift of flux tubes driven by the magnetostatic force 
(J X B per unit volume). More specifically, negative 
voltages have been observed by different groups^"— when 
recording the current-voltage characteristics at specific 
locations on the surface of the sample (central region). 
In addition, the striking effect takes place within a defi- 
nite interval of applied magnetic fields. 

Within such a complex scenario, it was recognized 
early on that the observations could only be under- 
stood if new dissipation mechanisms, additional to the 
flux flow phenomena were considered. In particular, a 
prominent role happens to be played by the flux-line 
cutting (crossing and recombination) between adjacent 
tilted vorticesi^ Nevertheless, certain facts still remain 
to be fully understood. Thus, the challenging problem of 
inhomogeneous electric fields, even changing sign along 
the specimen surface is hitherto open. On the other hand, 
important issues as the consideration of irreversible ef- 
fects related to the thresholds for flux depinning and 
cutting phenomena have not been reported yet. In this 
contribution, we investigate the influence of such mech- 



anisms on the establishment of critical negative current 
structures within the superconducting state. This is to 
be considered as a step forward for gaining knowledge on 
the processes that operate just previous to the dissipation 
regime. Within such a physical scenario, the concentra- 
tion of transport current towards the center of the sample 
and the appearance of negative flow at the surface will be 
predicted for a certain range of experimental conditions. 
To be speciflc, the application of different components of 
magnetic fleld, their sequence and characteristic values 
will be identifled as relevant issues for the observation of 
negative currents. 

The article is organized as follows. In Sec HI] we put 
forward the basic ideas about the theoretical approach 
used, that is a general critical state theory for type-II su- 
perconductors. An idealized slab geometry arrangement 
is proposed, aimed to introduce the lowest level of com- 
plexity for our purposes. Then, in Sec lIIII we perform a 
simplified analytical evaluation that allows to capture the 
main underlying physical matters. In Sec lIVI a quasi-3D 
statement of the problem is solved by numerical means. 
This is needed for the consideration of inhomogeneities 
over the sample. A discussion about the scope of our 
investigation within the problem of negative currents is 
given in SecfVl 



II. CRITICAL STATE APPROXIMATION 

A. Classical Maxwell equation approach 

The fundamental concept on which the critical state 
theory relies is that, in many cases, the experimental 
conditions allow to analyze the evolution of the system 
in the quasistationary regime. Thus, Ampere's law be- 
comes V X H = J, and determines the distribution of su- 
percurrents within the sample. When some external exci- 
tation (magnetic field and/or transport current) applies. 



the quasistationary evolution between successive equilib- 
rium states is ruled by Faraday's law V x E = — cJjB. 
Here the induced transient electric field is determined 
through an appropriate material relation J(E), and is 
used to update the profile of J. 

The material law encodes the mechanisms related to 
the breakdown of magnetostatic equilibrium, as well as 
the dissipation modes operating in the transient from 
one state to the other. In this sense, there have been 
a number of theoretical proposals, and among them (see 
Refsl^ andKLOl and the citations therein), here we choose 
the so-called double critical state model (DCSM) iii This 
approach allows a straightforward connection between 
the mesoscopic flux depinning and cutting phenomena, 
and the field equations for the coarse grained quantity 
J. On the one hand, the model establishes the criti- 
cal conditions |Jj_| < Jc± and |J||| < J^i that relate 



to (i) the maximum pinning force on the vortex lines 
(I J X B| = J±B < i^p.max) and (ii) to the maximum vari- 
ation of the tilt angle between vortices (notice that from 
Eq.Q one has da ex J\\)-— On the other hand, as it was 
thoroughly discussed in RefllOl the model also provides 
a rule that fixes the trajectory of the system through the 
dissipation excursion towards the new equilibrium state. 
Thus, corresponding to a very sharp transition from the 
superconducting state to some regime of high losses, one 
can argue that the final state current density verifies a 
maximum projection law relative to the transient electric 
field, i.e.: max (J • E) that ensures the fastest return to 
equilibrium. Notice that, in ID situations (infinite slab 
with a single component applied magnetic field) this is 
trivially verified because by symmetry one has J || E 
and both perpendicular to B. In other words, one has 
J± = sgn(Ej_)Jc±_ with E±_ standing for the component 
of E along the direction B x (J x B). In more general 
configurations, the maximum projection condition is not 
so simple due to the vectorial character of the problem. 
Thus, within the DCSM framework, the current density 
must transit from one state to another that fulfill the 
conditions \3 a_\ < Jci_ and |J||| < Jc\\. This may be ex- 
pressed by a relation of the kind J e A with A having 
a rectangular section of size 2 J^i x 2Jc± in this case. In 
summary, the critical state model in general 3D systems 
is posed by the system of equations 



B. Variational statement of the critical state 
problem 

From the mathematical point of view, the above prob- 
lem is equivalent to the incremental minimization of the 
functional (field Lagrangian) 
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- p • (V X H - J) 



d^ 



r. 



(2) 



Here, one introduces the variable p as a Lagrange mul- 
tiplier for enforcing Ampere's law. AH, on the other 
hand, represents the magnetic field increment for the 
time step under consideration. Additionally, the alge- 
braic condition J G A should be fulfilled in the mini- 
mization. 

Notice that the Euler-Lagrange equations for this vari- 
ational problem are 



dpC^O- 



V X H- J 



duC - d' 



dC 



0^ AioAH = -V X p (3) 



that identify p ~ EA<. From the mathematical point 
of view, here J is no longer a variable, but plays the 
role of a parameter to be adjusted in a direct algebraic 
minimization, i.e.: 



max (E • J) <=> max (p • J) 



(4) 



As one can see, both Ampere and Faraday's laws are 
included in the variational formulation, as well as any do- 
main A for the critical current material law. Technically, 
we emphasize that the shorter is the path step, the better 
agreement with the standard Maxwell equation formula- 
tion. On the other hand, numerical methods for discrete 
constrained minimization can be used as an alternative 
to the integration methods for Maxwell's equations, that 
happends to be very convenient for dealing with 3D prob- 
lems. 



C. Application: 2D and 3D slab geometry. 



V X E = -dtB 
V-B = 



V X H = J (B = ^oH) 

max J • E with J e A . (1) 



Notice that, as equilibrium magnetization is usually ne- 
glected in the critical state regime, B = /igH is used. 

The integration of the above system of equations sup- 
plemented by appropriate boundary conditions may be 
cumbersome, so that an alternative formulation has been 
frequently used, that is fully equivalent, and states the 
problem in a variational form. It is briefly explained in 
the forthcoming paragraphs. 



A further advantage of the methodology introduced 
above is that the constraint relation J g A allows a pic- 
torial representation that provides a useful tool for the 
understanding of the current flow structures that arise in 
the longitudinal configurations. Notice, that in our case 
(DCSM conditions) A may be depicted by the cylindrical 
region in FiglT] 

Regarding the specific details about the mathematical 
technique for obtaining the numerical solution of Eq. ^ , 
the interested reader is addressed to our RefllOl. There, 
we analyzed a number of situations that are easily trans- 
lated to the study within this work. Here, we will just 
mention that our proposal consists of transforming the 
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FIG. 1. (Color online) Better resolution in original paper. 
Panel (a): magnetic process considered in this work. A mag- 
netic field Hzo is applied perpendicular to the surface of a 
superconducting slab, that is later subjected to a transport 
current along the j/-axis and to an increasing field Hyo- Panel 
(b): the critical current restriction is represented by a cylin- 
drical region A around the local magnetic field axis (length 
2Jt,|| and radius Jc±)- Ap is the projection onto the slab (xy)- 
plane. 7 is the angle between the field and the z-axis. The 
lower inset depicts an element of the slab and the current 
backflow. 



volume integral over the whole space in ([2]), into a dou- 
ble integration over the sample's volume. Then, upon 
discretization, one solves for the distribution of current 
in a proper set of circuits, under the corresponding con- 
straints for the components of J parallel and perpendic- 
ular to the local magnetic field. Thus, when one chooses 
the infinite slab geometry depicted in FigHJ such cir- 
cuits are naturally defined by a collection of current layers 
within the xy-plane, and each carrying a current density 
given by [Jx(zi), Jy(zi)], with z.j the position of the layer. 
Notice that, owing to the planar translational symmetry, 
a z-component of J may be ruled out. 

We emphasize that the slab geometry enables to study 
the appearance of the focused physical phenomena with 
the least mathematical complication, i.e:. this config- 
uration allows to clearly determine the mechanisms re- 
lated to the negative currents. As a main fact it will 
be established that, when building the parallel configu- 
ration, the response of the superconductor depends on 
the limitations for the current density established by the 
depinning threshold Jc±, on the orientation of the local 
magnetic field, and eventually on the flux cutting restric- 
tion Jc\\ . This is easily understood at a qualitative level 
just by glancing at FigHJ The critical current restric- 
tion is given by the region Ap that is the intersection 
between the cylinder A and the ccy-plane, where the cur- 
rent flows. For moderate values of the angle 7 between 
the local magnetic field and the z-axis, Ap is an ellipse 
of semi-axes Jr\ and JK with 






Jc±/cosj = J,^^Hl + Hl + Hl/H, 



(5) 



An increase of the in-plane magnetic field component 
will result in a tilt of the cylinder, by an increase of 
the angle 7. Note in particular that, initially, the max- 
imum value of the in-plane parallel current density, JK , 
grows with the angle 7, independent of Jell (which is, 
thus, absent from the theory) until the maximum value 

is reached. 



Jl^ 



J", 



Then, the ellipse is truncated 

and eventually would be practically a rectangle of size 
2 Jc|| X 2Jci_ when 7 -^ 7r/2. Outstandingly, for large val- 
ues of X = Jc\\/ Jci. (long cylinders), the critical current 
along the parallel axis, J^y increases more and more as 
the weight of Hzo decreases and, furthermore, this quan- 
tity is always beyond the individual values Jc_l and J^n 



III. SIMPLIFIED ANALYTICAL MODEL 

Here, we show that some of the experimental features 
that will be obtained later on, from numerical calcula- 
tions, may be already predicted by a simplified analytical 
model. Let us consider the excitation process depicted in 
Fig H] for the particular case H^o = 0^7 = 7r/2 (the re- 
gion Ap is a rectangle with axes defined by the directions 
parallel and perpendicular to H). 

A. Governing equations 

Ampere's law takes the following form for the infinite 
slab geometry considered in this work 



dHy 

dz 



Jx 



dz 



Jy 



(6) 



On the other hand, following the theory issued in Ref. Ilfl. 
one can show that such expressions may be transformed 
into the polar form 



^ dz^ "^11 



dH_ 

dz ^-^^ 



(7) 



with H ~ yH'x + H^ the modulus of the magnetic field 

vector, and a = a.ta.n{Hy/Hx) the angle between such 
vector and the x-axis. 

Now, the thresholds of flux depinning and cutting im- 
ply the in-plane conditions 



|J|]'l<J,^i(7 = ^/2) = Je|| ; \Jl\<J, 



c± ■ 



(8) 



It is apparent that, in general, Eq.([7]) and the condi- 
tions in Eq.® would not straightforwardly lead to the 
solution of the problem. Typically, one should also use 
Faraday's law, either by explicit introduction of the re- 
lated electric fields (as in Ref. Iffh . or by our variational 



statement. Nevertheless, in this case (7 = 7r/2 => J^, = 

Jc||), the resolution noticeably simplifies. In fact, for the 
situation considered, we will have a combination of the 
cases Jn = 0, ± J^ii and J^ = 0, Jc± and integration of 
Eq.(l7]) is straightforward. For further mathematical ease, 
we will also consider Jc\\ and Jc± to be field independent 
constants in this work. 

The following normalization, based on the physical pa- 
rameters that define the problem, will be used; j = 
3/Jc±, h = H/Jcj_a and z = z/a (a is the thickness 
of the slab). The origin of coordinates will be taken at 
the center of the sample. 

Following the notation introduced in Rcf. |ll| we will 
refer to different zones within the sample that are, in 
brief, macroscopic regions where well defined dissipation 
mechanisms occur. Inserting our normalized units, there 
can be T zones, where only flux dcpinning (transport) 
occurs (j|| ^ ,j± = ±1), C zones, where only flux cut- 
ting occurs (j|| = ±XiJ± — 0), and CT zones where 
both transport and cutting occur (j|| = ±x,j± = il)- 
Finally, one will have O zones where neither flux trans- 
port nor cutting take place (j|[ = , jj^ = 0). Introducing 
these possibilities in Eqs.© and ([5]) one gets the follow- 
ing cases for the incremental behavior of the magnetic 
fleld in polar components 



dh = 



(0,C) 

±dz (T,CT) 



da ~ 







(0,T) 



±ix/h)dz (C,CT) 



(9) 



and all that remains for obtaining the penetration profiles 
is to solve successively (integrate) for h and a with the 
corresponding boundary conditions (evolutionary surface 
values ho,ao). The case selection has to be made ac- 
cording to Lenz's law. We note in passing that further 
specification related to the sign is usually included in the 
notation. Thus, a T+ zone will exactly mean dh ~ +dz. 



B. Magnetic process 

1. Application of current 

To start with, the application of the transport current 
along the y-axis produces a T-)_ zone 



dh = dz ; da = 

4 



(10) 



hxo 



1 



0, 



that penetrates from the surface until the point where h 
equals 0, i.e.: Zpo = 1 — Itr- In our units, Zpo = 0.5 for 
Itr ~ hxo ~ 0.5. An O zone appears in the inner region 
< z < Zpo as far as Itr < 1- 



2. Application of parallel field: initial steps 

The above situation remains valid until hyo is applied. 
Then, upon increasing hyo, flux line rotation starts on 
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FIG. 2. (Color online) Penetration of the magnetic field com- 
ponents and rotation angle in the longitudinal transport ex- 
periment {H^o = 0) for a superconducting slab of thickness 
2a, as calculated from Eq.((7]). The zone structure induced by 
increasing the field Hyo is marked upon some of the curves. 
The dashed line corresponds to the unstable regime (see text). 
Dimensionless units for h and z are defined in the text. 



the surface and the perturbation propagates towards the 
center in the form of a C_T+ zone defined by 

dh = dz ; da = — x/^ dz 

4 (11) 

h = hxo + z — 1; a = ao + xln[l + (z — l)/^o] , 

that covers the range z^ < z < 1, defined by a = => 
z^ = 1 -|- ho[exjp{—ao/x) — 1]- The former T+ zone is 
pushed towards the center and occupies the interval z^ < 
z < z; 



with 



core < z < z,. 



= 1 — /iQ. Finally, an O zone fills the 



The upper panes of Fig[2] sketch the above described 
structure marked on the cartesian components of the 
magnetic field vector. The transition points between the 
different regimes are highlighted. 



3. Application of parallel field: instability at the center 



The 0/T+/C_T+ structure remains valid until the 
condition z^ = <=> /iq = 1 is reached, i.e.: the modulus 
of h penetrates completely within the interval < z < 1. 
Then, the O zone disappears, and a T+/C_T+ struc- 
ture fills the sample. We emphasize that this configura- 
tion becomes unstable owing to the boundary condition 



hx{z = 0) =0 that is dictated by the symmetry of jy 
around the center. Thus, corresponding to the even be- 
havior oi jy(z), hx{z) must be an odd function. In phys- 
ical terms, flux vortices penetrate from the surface with 
some orientation given by the components of the vector 
{hx,hy). Owing to the critical condition for the pene- 
tration of the field dh/dz = 1, as soon as the modulus 
reaches the centre, flux rotation must take place there. 
This is needed for accommodating the vector to the con- 
dition h(z = 0) = (0,/iy(z = 0)) = (0,/i(z = 0)). On 
the other hand, as the angle variation is determined by 
the value of J^n, a jump is induced at the centre, i.e.: 
a(z = 0) — > 7r/2, and the related instability may be vi- 
sualized by a critical C+T+/T+/C_T+ profile (dashed 
line in Fig. [2]) in which the field angle decreases from its 
surface value ao to in the C_T+ region, then keeps null 
within the T_|_ zone, and suddenly increases to the value 
7r/2 in the inner C+T^ band defined by 



dh = dz 



x/hdz 



hxo 



da 

4 (12) 

1; a = V2-xln[l+z/(/io-l)]. 



In fact, a C+T+/C_T+ structure is stabilized 
with the intersection between regions at the point 



a 



+,+ 



(Zv) 



a '+(zv)] given by 



1 



ho 



y/ho{ho — l)exp[(7r/2 — ao)/x\- Note that, upon further 
increasing hyo the point Zv follows the rule Zv(/ij,o ~^ 
cxd) — >■ {1+hxo/x)/'^- All these features have been marked 
in the lower pane of Fig. [2] . 



4- Physical phenomena 

The previous results allow to identify the following 
properties as hyQ is increased: (i) the appearance of a sur- 
face layer with negative transport current density (mind 
the slope of hx in Fig. ([2]) in view of Eq.®), and (n) the 
applied magnetic field re-entry as related to the inner 
C+T+ zone. These features will be confirmed along the 
forthcoming paragraphs, where the numerical solution of 
the problem is presented. Additionally, we will show that 
the inclusion of a third component of the magnetic field 
[hzo 7^ in what follows) allows to unveil further details 
reported in the literature. In particular, the occurrence 
of the negative current phenomenon at specific locations 
on the surface of the sample and for a given range of 
applied magnetic field will be understood within the 3 
dimensional scenario. 



IV. NUMERICAL RESULTS 

Next, we detail the results obtained numerically for 
different material laws, as related to the selection of the 
critical current region A. We restrict the plots to the lim- 
iting cases jell ~^ oo and j^n = 1 (infinite and unit aspect 
ratio, or also named after T and CT states for obvious 



reasons). The information of interest for intermediate 
values is straightforwardly interpolated. 

From the technical side, we must point out that when 
minimizing L [see Eq.([2])] a slightly smoothed version of 
the cylindrical region A has been considered by means of 
a superelliptic relationi^ given by 



ji" + (jii/x)'"<i 



(13) 



with n = 4. This statement performs with a high sta- 
bility from the numerical point of view. Figs. [3] and |3] 
display the main features obtained for the longitudinal 
transport experiment when the third component of the 
magnetic field (hzo) is incorporated. First, we will ana- 
lyze the properties of the field (h(z)) and current density 
(/(z)) profiles, for a longitudinal configuration built in 
the fashion described in Fig.([T|). Fig. ([3]) contains the be- 
havior of these quantities as hyo (applied parallel field) 
is increased, subsequent to the application of the trans- 
port current. This is done for a low and a high value of 
the perpendicular magnetic field hzo- For the meaning 
of low and high recall that, along this work, the units 
arc relative to the characteristic penetration field value 
H* = Jc±a 



A. Field and current density penetration profiles 

The curves in Fig. ([3]) provide a basic mapping of the 
physical conditions in which negative currents occur. For 
completeness and for quantitative purposes, we have in- 
cluded both the field and current density profiles. Nev- 
ertheless, they are related by Ampere's law (Eq.®) as 
one can easily check at qualitative level, i.e.: in the slab 
geometry J^^y is the slope of Hy^^ respectively. 

Recall that negative values for the transport current 
density jy are neither obtained for the T or CT states 
when hzo is high (hzo > 50) until extreme values of the 
longitudinal field {hyo > 1000) are reached. On the con- 
trary, one can early find negative current flow for both 
cases when h^o = 1. If j\\ is unbounded (T states) the 
jy{z) structure becomes rather inhomogeneous as hyQ in- 
creases and takes the form of a highly positive layer in the 
center shielded by a prominent negative region. When jy 
is bounded (CT states) one observes a negative layer at 
the surface that eventually disappears when hyQ increases 
more and more {hyQ > 50). 

Some fine structure details are also to be noticed: 
(i) for the magnetic process under consideration, the 
partial penetration regime in which the fiux free core 
progressively shrinks to zero (curves labelled hyQ = 
0.005 • • • 0.845) is practically independent of the criti- 
cal current model (region) used, (ii) the peaked struc- 
ture oi jy{z) for the T-states at h^o ~ 1 (curves labelled 
hyo = 10 • • • 300) is accompanied by a similar behav- 
ior in jxiz) that relates to a subtle magnetic field reen- 
try phenomenon in hy{z) as outlined in the plot, (iii) 
the negative values oi jy{z) are obtained for smaller and 




FIG. 3. (Color online) Profiles of the magnetic field components hx[z,hy{a)] and hy[z,hy{a)], including their corresponding 
current-density profiles jy[z,hy{a)] and jx[z,hy{a)] for the T-state model (left) and the superelliptic-CT model (right). The 
magnetic dynamics for high (top) and low (bottom) perpendicular magnetic field hzo is shown. Different scales are used to 
visualize the intricate behaviour and to emphasize the appearance of negative currents at the surface. The curves are labelled 
according to the longitudinal magnetic field component at the surface of the slab, hy{a) = hyo. 



smaller hyo as h^o also decreases. In fact, negative val- 
ues can happen even for the partial penetration regime 



{h. 



yO 



< 



0.845) when hzo tends to 0, in accordance with 



the analytical model presented before (Sec. (|III|) ). 



B. Experimental quantities 

For a closer connection with real experiments, we have 
also calculated the sample's magnetic moment M as a 
function of the longitudinal field. The fingerprints of neg- 
ative current flow will be identified. Fig. |4] displays the 
magnetization process of the slab as a function of the 
applied longitudinal field amplitude hyQ. M (in units of 
Jc±ol^), as well as the transport current density at the 
center jy{0) and at the surface jy{a) are displayed. 

Several features are to be identified: 

(i) Unhmited growth of Mx{hyo) and jj;(0) occurs for 
the T-states, in which j|| is unbounded. On the other 
hand, the appearance of a peak structure in M^ihyo) 
correlates with a maximum value of the transport cur- 
rent density at the center of the slab for the CT states. 
The obtained maximum value jy''°'^'{Q) — 1.2968 corre- 
sponds to the optimal orientation of the region A in 
which the biggest distance within the superelliptic hy- 
pothesis is reached. Such situation is sketched in FiglS] 
and one may check the numeric result from the expres- 
sion 



maxj^|| = (l + x'"/("-'^)^"-')/'" 



(14) 



when the choices x = 1 j '* = 4 are used. The above for- 
mula is obtained from Eg. p^ after straightforward cal- 
culations. Notice that, as a limiting case, it produces the 
expected value 2^^^ for the diagonal of a perfect square 
(i.e.: rn- oo in Eq. (fT3)) ). 

(ii) For the unbounded case, in the low hzo regime, the 
negative current density at the surfaces stabilizes towards 
the value jy{a) = —0.5. 

(iii) For the bounded case, and moderate or low h^o, 
the transport current at the surface stabilizes towards the 
value jy{a) = 0.422 whether or not it has been negative 
along the ramp of applied longitudinal field. 

(iv) As a general rule, the smaller the value of h^o, the 
sooner the negative transport current is found. In the CT 
case, this also increases the range of longitudinal field for 
which negative values are observed. 



V. DISCUSSION 

Within the previous sections we have displayed a num- 
ber of cases in which negative transport layers are pre- 
dicted for type-II superconductors if the longitudinal field 
configuration (J || B) is induced by some external pro- 
cess. Our theoretical investigations allow to identify the 
following relevant aspects for the appearance of such phe- 
nomenon: 

(i) the physical mechanisms of flux cutting and depin- 
ning, that may be modelled by the thresholds for the 
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FIG. 4. (Color online) Top: components of the magnetic moment of the slab as a function of the applied longitudinal magnetic 
field hyo — hy{a) for the T-state model (left) and the superelliptic-CT model (right) with different values of the perpendicular 
magnetic field component h^o as labelled by arrows on each curve. The magnetization dynamics at low fields hyo is shown in 
the 3D-insets. Notice the different scales. Bottom: the corresponding evolution of the transport current jy at (z = 0) and 

{z = a). 




xy - plane 



FIG. 5. (Color online) Side view of the critical current region 
A rotated by an angle 7 that produces a maximal parallel cur- 
rent at the xy-p\a,ne (this plot is a specific longitudinal section 
of Fig[T]). The precise orientation takes place for a definite 
value of the applied magnetic field hyo- The superelliptical 
region considered in this work has been plotted together with 
the strictly CT model. 



components of J parallel and perpendicular to the local 
magnetic field. In this sense, we have shown that neg- 
ative currents are much easier observed for materials in 
which Jc|| and Jc± are independent and J^-u ^ Jc±. On 
the other hand, additional calculations (not displayed) 
imply that alternative ansatzs as the isotropic model 



(J 



Jl 



< Jt, i.e.: the difference between the mech- 



anisms responsible for the thresholds J^n ^-nd Jc± are 
not relevant and the region A is a circle) cannot predict 
such behavior. However, if some interaction is allowed 
between the cutting and depinning thresholds [moder- 
ately smoothed Jj^(J||) relation, i.e.: the region A is a 
superellipsoid] , the negative current flow will occur for 
some range of conditions. 

(ii) The sequence in which the external excitation com- 
ponents are applied to the superconductor. Thus, as 
one could expect from the idea that in the critical state 
all changes proceed from the surface toward the centre 
of the sample, negative current structures are enhanced 
when the applied magnetic field is applied after origi- 
nally building a subcritical transport profile. Recall that 
in the situations depicted in Fig|3l the initial current flow 
is compressed more and more, even until compensating 
negative values are needed at the surface, for maintain- 
ing the biased transport current. Along these lines, we 
should comment that when simulating experiments in 
which the transport current is applied subsequent to the 
field, our theory does not predict negative fiow values at 




FIG. 6. (Color online) Penetration of a magnetic field parallel 
to the axis of a finite superconducting cylinder (half longitu- 
dinal section is shown for symmetry reasons). The compo- 
nent of the magnetic field perpendicular to the lateral surface 
is visualized by a set of arrows with normalized lengths. The 
dashed line represents the symmetry axis. 



all. On the contrary, in such cases, what one gets is a 
compression of the original field penetration profile, until 
the increasing transport current leads to dissipation. 

Additional physical considerations can be done so as 
to cover the full experimental scenario. In particular, 
though our analysis has been done within the infinite 
slab geometry, one can straightforwardly argue about the 
extrapolation to real experiments. Thus, the inclusion of 
the third component of the magnetic field H^o relates to 
the last relevant aspect: 

(iii) The importance of the finite size effects. Notice 
that from our numerical calculations, one can predict 
that negative currents should be more prominent in those 
regions of the sample where the component of H perpen- 
dicular to the current layers is less important. This will 



occur at the central region of the sample's lateral sur- 
face, where end effects are minimal. Thus, considering 
that a real sample in a longitudinal configuration will be 
typically a rod with field and transport along the axis, 
the above idea is straightforwardly shown by plotting the 
penetration of an axial field in a finite cylinder. This has 
been done in Fig|6l The plot shows the distortion of the 
magnetic field, shielded by the induced supercurrents. 
Just for visual purposes, we have superimposed the hori- 
zontal component of the magnetic field, along the lateral 
side of the cylinder. It is apparent that the normal com- 
ponent of H will be enhanced close to the bases and tend 
to zero at the central region. Thus, inhomogeneous sur- 
face current densities, with negative flow at the mid part 
should be expected. 



VI. CONCLUDING REMARKS 

In this article, we have shown that the counterintu- 
itive effect of negative current flow in type-II supercon- 
ductors may be predicted and quantifled by means of the 
critical state theory. For restricted situations (infinite 
slab geometry and fields parallel to the surface) , the pre- 
diction may even be done within a simplified analytical 
model. Three dimensional effects may be incorporated 
by numerical methods when a third component of the 
magnetic field, perpendicular to the surface of the slab 
is considered. The analysis of this situation has allowed 
to conclude that negative transport is enhanced for su- 
perconductors in which the flux cutting barrier is much 
above the depinning value, and at those regions of the 
sample where the magnetic fleld is basically parallel to 
the surface (central part in cylinder geometry) . 

Extensions of this work are planned along two lines: 
(i) the actual evaluation of the longitudinal problem in 
flnite length samples (fleld and transport along the axis 
of a rod) , and (ii) the extrapolation of the current den- 
sity profiles beyond the critical state threshold, so as to 
include the dissipation effects. 
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Recall that the relation {x/a) " + (j//fe) " = 1 represents a 
family of closed curves (superellipses) with semi-axes a and 
h, ranging from an ellipse to a rectangle as n follows the 
sequence n = 1, 2, 3, . . . cxd . This was exploited in RefflOl 
for comparing different critical state material laws. 



